Degradation versus self-assembly of block copolymer micelles 
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The stability of micelles self- assembled from block copolymers can be altered by the degradation of 
the blocks. Slow degradation shifts the equilibrium size distribution of block copolymer micelles and 
change their properties. Quasi-equilibrium scaling theory shows that the degradation of hydrophobic 
blocks in the core of micelles destabilize the micelles reducing their size, while the degradation of 
hydrophilic blocks forming coronas of micelles favors larger micelles and may, at certain conditions, 
induce the formation of micelles from individual chains. 
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I. INTRODUCTION 

Physicochemical properties of block copolymers of- 
ten determine their function in many useful applications 
ranging from biotechnology and drug delivery to paint- 
ing and oil extraction^-^. The possibility to control es- 
sential properties of block copolymers capable to adapt 
the behavior to changes in the environment is thus a 
challenging task. One of the main properties of am- 
phiphilic diblock copolymers in solution is their ability 
to self-assemble in micelles^ composed of a hydrophobic 
core surrounded by a hydrophilic coronaii. Such com- 
pact finite size aggregates can encapsulate hydropho- 
bic agents in their cores^. In particular, this loading 
capacity of block copolymers can be used for selective 
transport of hydrophobic nanoparticles and lipophilic ac- 
tive molecules to specific targets and through the cell 
membrane^. 

However, the use of block copolymer micelles for tar- 
geted drug delivery implies also the necessity to control 
the release of the transported particles from the cores of 
micelles, for example by external stimuli. In turn, the re- 
lease process is closely related to thermodynamic stabil- 
ity of micelles. Micelles assembled from block copolymers 
can be relatively stable^. This hinders the release of ac- 
tive components^ and thus limits their use for biomedical 
applications. 

Degradable polymers^, provide for additional degree of 
freedom allowing to control the longevity and stability 
of block copolymer micelles. Degradation of the polymer 
backbone may change significantly the thermodynamics 
of block copolymer self-assembly and thus stability of the 
micelles. Tuning the rate of degradation would allow for 
modulation of the thermodynamic stability of micelles in 
large extent. In addition, tumor tissues have tendency to 
selectively accumulate polymers. This effect is known as 
the enhanced permeation and retention (EPR)iS and is 
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attributed to larger size of the pores in blood vessels of 
tumor tissues. Polymers can interact with the cell mem- 
branes and reticuloendothelial system (RES)^^, what can 
potentially increase their cytotoxicity^^. From this per- 
spective, the degradation of polymers up to metabolites^^ 
may solve biocompatibility issues and excessive accumu- 
lation in tissues. 

Usually biomedical applications require long- 
circulating delivery vectors with constant release 
for days^^. This implies that polymer degradation in 
such systems is much slower than the time required 
to reach the thermodynamic equilibrium. Kinetics of 
micelle formation from monomers can be rather fast. 
Characteristic times of exchange of oligomers between 
micelles and the bulk and the relaxation do not exceed 
milliseconds^^, while for longer block copolymers, the 
characteristic time is of order of minute&lSJi. Thus, in 
such systems the self-assembly in micelles is a quasi- 
equilibrium process with a fixed length distribution of 
the blocks, that changes with time. Micelles of degrad- 
able block copolymers can assemble and reassemble 
changing their internal structure and shape^i^^ thus 
changing the loading capacity of the cores. The balance 
between steric repulsion of hydrophilic blocks in the 
coronas of micelles and interaction of hydrophobic 
blocks forming part of the cores defines the finite size 
of micelles. Degradation of hydrophilic blocks^i2£ would 
then lead to destabilization of micelles leading to fusion 
and formation of larger cores. In turn, degradation 
of hydrophobic blocks in the cores^i"— may induce 
splitting the cores of micelles. 

Using a scaling theory we study the interplay between 
the degradation of one of the blocks of block copoly- 
mers and equilibrium self-assembly and re-assembly of 
block copolymers into micelles in case of both corona and 
core degradation. We present a scaling theory of block- 
copolymer self-assembly into micelle o^^'^^ , and discuss 
the degradation kinetics of the polymer blocks. How- 
ever, we note that scaling theories have limitations, (i) 
Scaling theories in principle only apply to (infinitely) 
long chains, (ii) the theory does not include kinetics; 
(iii) it does not include morphological changes. We 
assume two types of degradation mechanisms, (i) ran- 
dom chain scission mechanisnt^, which implies stochas- 
tic process of a polymer chain division at random posi- 
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FIG. 1. Quasi-equilibrium degradation of micelles. 

dation induces growth of micelles. 



Core degradation induces the disassembly of micelles while corona degra- 



tion. This mechanism causes a gradual decrease in the 
number average molecular weight. Such random divi- 
sion of a chain is usually observed in hydrolytic degra- 
dation of polyesters"^^'^^ or certain polyamides^^^. (ii) 
End-evaporation mechanis m^^i^^ , when the monomers 
are gradually detached from the ends of the polymer 
chain. This mechanism is typical for enzymatic degra- 
dation. The degradation of the blocks shifts the critical 
micellar concentration (CMC) and changes the size dis- 
tributions and the average sizes of the micelles (Figure 
[1]) . The degradation of the chains can induce or suppress 
the self-assembly process. 



II. SELF-ASSEMBLY 

Diblock copolymers composed of soluble and insolu- 
ble blocks can spontaneously self-assemble into micelles^. 
The micellization is the entropy driven process, the en- 
tropy of mixture of individual chains in the solution 
is balanced by the tendency of insoluble blocks to re- 
duce contacts with the solvent in the cores of micelles. 
Thus, the stability of the micelles is defined by the 
energy of insoluble blocks forming compact cores and 
the steric repulsion of soluble blocks in the coronas of 
micelle a^^'^^ . Characteristic equilibration time of block 
copolymers depends exponentially on the length of the 
insoluble bloc k^^'^^ . In the cases of long insoluble blocks 
and high interfacial tensions^l the relaxation time may 
be hours, and thus non-equilibrium kinetics should be 
considered. However, we focus on situations when the 
hydrophobic block is relatively short or the interfacial 
tension is low. In this case the relaxation time is small 
(milliseconds)^^ and thus the kinetics of self-assembly is 
much faster than the degradation. If the degradation 
process is slow enough to allow for the equilibrium assem- 
bly of block copolymers into micelles, the scaling model 
of equilibrium self-assembly can be applied. We assume 
that only one block, either insoluble or soluble, is degrad- 



able, so that the total number of block copolymers does 
not change with time. In addition, we consider spherical 
micelles, while morphological transitions are not consid- 
ered. We denote Cp the number density of aggregates 
comprising of p copolymers, where p is the aggregation 
number. The total free energy per unit volume of the 
solution of copolymers at a given time is 



kT 



c{n) ln(^A3)dn 
e 

(1) 

where the first term is the entropy and Fp is the free 
energy of a micelle comprising p copolymers, A is the 
de Broglie wavelength, p — I corresponds to individual 
block copolymers contributions, while the entropy of free 
fragments is taken into account in the last term. Degra- 
dation of blocks provokes the detachment of fragments of 
different lengths floating in the solution and the last term 
takes into account the entropy of fragments where c(n) is 
their length distribution function. The total number of 
monomers in the self-assembly and the degradation pro- 
cess is conserved. This is reflected in the conservation of 
mass condition, 



^pcp = LP 



(2) 



where ip is the total copolymer concentration. Since we 
consider the degradation of one block, the total number 
of copolymers in the solution is not changed with time 
and the degradation does not affect this condition. Min- 
imization of the free energy ([T]) subject to the constraint 
gives the quasi-equilibrium distribution of the copoly- 
mers in the micelles^. 
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This expression describes the distribution of micehes of 
degradable copolymers at each time. 

Exphcit form of the free energy Fp is defined by the 
molecular structure of copolymers forming a micelle and 
is the sum of the corona and the core contributions, 



Fp = F™™"'' + (4) 

The free energy of individual chains, p = 1, is also de- 
scribed by this expression, where the corona term trans- 
forms in the entropy contribution of a linear chain and 
the insoluble block gives the corresponding core contribu- 
tion. The exact expressions of Fp°^°^^°- and Fp"^*^ depend 
on degradation mechanism and are functions of time. 

In the following, we denote the chain length of the 
soluble block as N and the chain length of the insoluble 
block as Nc- We assume sufficiently long soluble chains, 
TV ^ 1, and the monomers of both blocks being of the 
same size a. 



III. DEGRADATION KINETICS 

We consider random chain scission and end evapora- 
tion mechanisms and degradation of soluble and insoluble 
blocks. 



1. a) Core degradation 

Insoluble blocks tend to avoid contacts with water and 
favor assembly into the core of the micelle. The degrada- 
tion of insoluble blocks makes block copolymers more sol- 
uble and this would shift the equilibrium towards smaller 
micelles. Assuming that the core of a micelle is a dense 
and homogeneous sphere formed by p copolymers, the 
free energy of the core is given by-^ 

o 2 ^ 

Fp™''- = kT{36TT) i aNipi + (6) 

807Vc' 

where a is the surface tension of a sphere of radius 
Rc = (3/(47r)pA^c)^^"^a and the second term describes the 
elastic contributions arising from the stretching of insolu- 
ble blocks in the coreSE. The " effective" surface tension of 
the core a describes implicitly the fact that the insoluble 
chains tend to avoid contact with the solvent by forming 
a dense core. The lengths of the insoluble blocks Nc de- 
pend on time and for each time t the length distribution 
is given by equation of random scission degradation ([5]). 

Steric repulsion in the corona of the micelles, formed by 
soluble blocks, penalizes the formation of large micelles. 
Since the soluble block does not degrade, the corona of 
micelles is composed of chains of equal length N . The 
partition function of a monodisperse star Zp yields in the 
form^^ 



A. 1) Random chain scission 

Random scission mechanism implies homogeneous dis- 
tribution of splitting points along the chain. At a given 
time a chain is divided into random parts. Thus, the dis- 
tribution of fragments of different lengths P(n, t) is given 
by^o 



Zp oc iV>" 



(7) 



where N is the length of the arrc^i. 7p are the universal 
critical exponents of the star polymers^. The numerical 
values of 7p are known exactly for a wide range of 
and in the range < p < 200 can be interpolated by the 
power law expression 7p = 1 — 0.0893(p — 1.5)^'^^. With 
this, the free energy of polydisperse corona is given by 



nP{n,t)+ / dyP{y,t) (5) 



1 dP{n,t) 
K dt 



The negative term refers to the loss of chains of length 
n, and the positive term describes the gain of chains of 
length n due to the degradation of longer chains (with 
lengths more than n) at a given time t. The coeffi- 
cient K takes into account the fraction of chains that 
degrade simultaneously at each time, practically defin- 
ing the timescale. Initial distribution of chains at i = 
is assumed to be monodisperse. This equation is solved 
numerically by considering the integral in the right hand 
side Q{n,t) = dyP(y,t) a numerical function. This 
function is calculated for a given time as a function of n 
and then the equation ([5|) is solved numerically. Typi- 
cal length distribution functions are shown in Figure [21 
left column. Initial homogeneous distribution gradually 
disperses and shifts to small n. 
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-kTliiZp = -kT{jp-l)\iiN (8) 



Eqs. ([6]) and ([8]) define the free energy of micelles of p 
copolymers and the free energy of individual chains for 
p = 1. It allows to calculate the size distribution of mi- 
celles as a function of time The results are present 
in Figure [3^) The degradation of the core starts when 
the micelles are formed (concentration above CMC). The 
micelles gradually decrease in size and disassemble. The 
degradation rate k is related to the time step in the degra- 
dation equation ([5]). The chosen value k = 0.0003 is low 
enough to insure gradual changes in the size of the mi- 
celles, if the rate is higher, the micelles would disappear 
faster. In a real experimental situation this parameter 
connects the time step with real time, e.g. k = 0.0003 
signify that only three of ten thousand chains disassoci- 
ate at one time step. In experiments this parameter may 
vary in a wide range. 
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FIG. 2. Length distributions in the random scission (left column) and the end evaporation (right column) mechanisms for 
different degradation rates k, and initial length 200. 



2. b) Corona degradation 

Situation is different when the soluble blocks can de- 
grade while the insoluble blocks are stable. The cores 
of the micelles are formed by insoluble blocks and the 
core contribution has the same form as in previous case, 
eq. ([ni, but Nc remains constant. In turn, soluble 
blocks forming corona now can degrade with time and 
the corona contribution changes. Coronas of micelles are 
formed by polydisperse arms with the length distribution 
P(n, t) given by The partition function of a polydis- 
perse star Zp yields in the form— 



Zp oc nJ'-^--'nJ'-'~^--' . . . n^ri'^^n?-' (9) 

where ni < n2 < . . . < n^^i < rip are the lengths of the 
corresponding arms, sorted in ascending order^^. This 
partition function leads to the corresponding expression 
of the free energy of the corona. 



7p_2)lnn2 + . . . + (71 - 1) In rip] 

This expression defines, together with eq. the free 
energy of the micelles and individual chains in the solu- 
tion (|4]). In fact, the scaling expression of the corona con- 
tribution of a micelle with one arm, p^°^°^^ = —l^T{ji — 
1) In corresponds exactly to the scaling expression of a 
soluble chain in a solutior^. 

The resulting size distribution of micelles ([3]) is shown 
in Figure ^jp) and c) for different times. The effect of 
corona degradation is opposite to the case of core degra- 
dation. Figure shows the corona degradation kinetics 
for the same initial conditions as in Figure [3^) when the 
micelles are formed for conditions above CMC. Degra- 
dation of corona induces self-assembly of micelles from 
initially homogeneous solution of individual chains. In- 
dividual chains associate into micelles and the distribu- 
tion of micelles grows. Consequent degradation leads to 
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FIG. 3. Time evolution of the normalized size distribution 
pcp. a) core degradation above CMC, a — 0.81, k — 0.0003 b) 
corona degradation above CMC, a = 0.81, k — 0.01 c) corona 
degradation below CMC, a — 0.75, k — 0.01. The initial 
lengths of the blocks are: A'^ — 200, Nc = 20, copolymer 
concentration ci — 10~®. 



the shrinkage of the corona, thus larger micelles become 
more stable and small micelles unstable. The size dis- 
tribution of the micelles becomes broader with time due 
to increased polydispersity. Further degradation of hy- 
drophilic blocks would lead to the growth of micelles, 
morphological changes^S of the micelle shape and conse- 
quent bulk phase separation. 

Another important effect of corona degradation is the 
possibility to induce self-assembly of micelles from the 
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FIG. 4. Time dependence of the critical micelle concentration 
(CMC). Parameters are the same as in Figure [3^) and b). 

solution below CMC (Figure [St)) The degradation in- 
duces formation of micelles from initially homogeneous 
solution first for small aggregation numbers. Consequent 
degradation increases the number of micelles, the poly- 
dispersity and the average size. Thus, the degradation 
influence the CMC, which is time dependent. Figure 2] 
It decreases with time in case of corona degradation (red 
curve) and increases with time in case of core degradation 
(blue curve). 
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FIG. 5. Time evolution of the aggregation number of the 
maximum of the size distribution, PmUX. Parameters are the 
same as in[3l 

The effect of degradation on micellization process can 
be summarized in the plot showing the position of the 
maximum Pmax of the size distribution Cp , Figure [51 In 
case of core degradation (blue line) the maximum of the 
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distribution moves to small numbers (see Figure |3^)) un- 
til the micelles disappear completely (dashed line) , while 
in case of corona degradation (red curve) the size of the 
micelles increases (see Figure [5)d)) until the morphology 
changes or phase separation occur. If the degradation 
starts below CMC, the self-assembly into micelles starts 
at a given time which depends on the rate of degradation 



B. 2) End evaporation 

Random division of a chain is not the only degradation 
mechanism of polymer chains. Degradation in certain 
chemical reactions and enzymatic degradation may lead 
to gradual decrease of the chains from the ends (chain- 
end-activated degradation). The process of loosing the 
monomers from the end is described by the following pro- 
cess, P{n,t+1)-P{n,t) = K{P{n + l,t)-P{n,t)) which 
can be written in the integral form similar to ([5]) as 




P 



FIG. 6. End evaporation kinetics of core degradation in com- 
parison with random scission. Parameters are the same as in 
Figure [3^) except the degradation rate, k — 0.03. 



-^^^^ = -P{n,t)+ r 6{y-n~l)dyP{y,t) (11) 

where S is the Dirac function. Using this equation in- 
stead of dni, one can obtain the time dependence of the 
free energy of micelles for core and corona degradation. 
Figure [5] right column presents the kinetics end evapora- 
tion degradation. The starting and the final chain length 
distributions are very close for both types of degrada- 
tion. That is why the micelle distributions shown for core 
degradation in Figure IH] also coincide. The difference is 
only seen for intermediate times. However, end evapo- 
ration degradation is much slower than random scission, 
thus the rates of degradation differ in this example 100 
times. 

In conclusion, scaling theory of quasi-equilibrium mi- 
cellization coupled with the degradation of the blocks 
demonstrates that the degradation of hydrophilic blocks 
can induce self-assembly of copolymers into micelles and 
increase the size of the micelles, while the degradation of 
the hydrophobic blocks destabilize the micelles, reduce 
the equilibrium size of the micelles and can lead to com- 
plete disassociation of the micelles. It is valid for random 
scission mechanism assumed for degradation mechanics 
as well as for enzymatic (chain-end) scission mechanism. 
These findings may suggest the ways of controlled self- 
assembly and destabilization of micelles by degradation 
of the blocks. Our model do not account for morphologi- 
cal transitions and we plan to study them in future with 
a more detailed microscopic theory^. 
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